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ABSTRACT 

We report new, sensitive observations of two z ~ 3-3.5 far-infrared-luminous radio galaxies, 
6C 1909+72 and B3 J2330+3927, in the 12 CO J = 1 - transition with the Karl Jansky Very 
Large Array and at 100-500 /im using Herschel, alongside new and archival 12 CO J = 4 — 
3 observations from the Plateau de Bure Interferometer. We introduce a new colour-colour 
diagnostic plot to constrain the redshifts of several distant, dusty galaxies in our target fields. 
A bright SMG near 6C 1909+72 likely shares the same node or filament as the signpost active 
galactic nuclei (AGN), but it is not detected in 12 CO despite ^20,000 km s _1 of velocity 
coverage. Also in the 6C 1909+72 field, a large, red dust feature spanning «500kpc is aligned 
with the radio jet. We suggest several processes by which metal-rich material may have been 
transported, favouring a collimated outflow reminiscent of the jet-oriented metal enrichment 
seen in local cluster environments. Our interferometric imaging reveals a gas-rich companion 
to B3 J2330+3927; indeed, all bar one of the eight z > 2 radio galaxies (or companions) 
detected in 12 CO provide some evidence that starburst activity in radio-loud AGN at high 
redshift is driven by the interaction of two or more gas-rich systems in which a significant 
mass of stars has already formed, rather than via steady accretion of cold gas from the cosmic 
web. We find that the 12 CO brightness temperature ratios in radio-loud AGN host galaxies are 
significantly higher than those seen in similarly intense starbursts where AGN activity is less 
pronounced. Our most extreme example, where L' C04 _ 3 / L' coi _ > 2.7, provides evidence 
that significant energy is being deposited rapidly into the molecular gas via X-rays and/or 
mechanical ('quasar-mode') feedback from the AGN, leading to a high degree of turbulence 
globally and a low optical depth in 12 CO - feedback that may lead to the cessation of star 
formation on a timescale commensurate with that of the jet activity, <10Myr. 

Key words: galaxies: active — galaxies: high-redshift — galaxies: starburst — submillimetre 
— infrared: galaxies — radio lines: galaxies 



1 INTRODUCTION 

High-redshift radio galaxies (HzRGs) are typically identified via 
their ultra-steep-spectrum radio emission (a < — 1 where S v oc 
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v a ) in flux-limited s urveys (e.g. iTielens. Milev & Willis! 1 19791 : 
Rottge ring et alj [l 994i) . Despite the inevitable and extreme youth 
of the radio jets that draw our attention to these distant galax- 
ies, they are associated with the most massive stellar popula- 
tions of any known galaxy class - and presumably the most 
massive black holes and host galaxies - out to the highest 
redshifts (Best, Longair & Rottgerin gHl998t iBlundell & Rawlingsl 
ll999l : ISevmour et all2007l) ~ 

In the submillimetre (submm; rest-frame far-infrared, FIR) 
regime, HzRGs were first explored using relatively primitive 
submm detect ors such as UKT14 on the James Clerk Maxwell 
Telescope (e.g. Dunlop et al. 1993). Later, with the advent of sen- 
sitive submm cameras, their submm luminosities (and hence dust 
mass es) were found to be a strongly increasing function of red- 
shift dArchibald et al.ll200ll : lReuland et al.ll2004l) . even beyond the 
peak epoch of activi ty for submm galaxies (SMGs, z ~ 1-3 - 
IChapman etai]|2005l) . 

As rare, massive systems, HzRGs are often employed as 
signposts to what are expected - based on our understand- 
ing o f how cosmic structures form and evolve (e.g. iDavis et al] 
1985) - to be over-dense regions of the early Universe. Surveys 
found excesses of various galaxy types around Hz RGs, includ- 
ing Lyman-break galaxies, Lya emitters and SM Gs (Ivison et al. 



Papadopoul os et alj2005h or detections have been of a tentative na - 



200d: iMilev et alj|2004 iGreve et alj|200l IVenemans et al.ll2007l : 
Overzier et al.ll2008l : lHatch et alj|201 ll) . IStevens et all d2003l) pre- 



sented submm imaging of seven HzRGs, several of which appear 
to contain extended (~5-20 arcsec, ~35-150kpc) dust emission, 
co-spatial with similarly ex t ended UV emission in s everal cases 
l lHippelein & Meisenheimeri 1 19931 ; lHatch et alj 120081) . This sug- 
gests that obscured starbursts in these over-dense regions at z Z 3 
may differ from the compact ( <0.5 arcsec, or <4 kpc) events seen in 
local ultraluminous IR galaxies (ULIRGs) and, indeed, as seen by 

high-resolution 12 CO and radio continuum i maging of the general 2 OBSERVATIONS 

' et al.ll2008l : lBiggs & Ivisor] 



ture (e.g. Greve. Ivison & Papadopoulos 2004; Kl amer et al.l2005h . 

Here, we present new observations of two HzRGs, 
6C 1909+72 (als o known as 4C 72.26 and TX S J 1908+7220) and 
B3 J2330+3927 dPentericci et al.l2000LlDB03l) using JVLA, the In- 
stitut de Radioastronomie Millimetrique 's Plateau de Bure Inter- 
ferometer (IRAM PdBI) and HerschM dPilbratt et al.ll2010t). The 
PACS and SPIRE instruments jGriffin et al.ll2010l : |Poglitsch et al.1 
|2010|) aboard Herschel offer unprecedented sensitivity, with reso- 
lution well-matched to ground-based predecessors such as SCUBA 
dHolland et alJI 19991) . They cover the decade of wavelengths from 
70 to 500 with A/ A A ~ 3, which allows us to probe the peak 
of the spectral energy distribution (SED) of a dusty galaxy out to 
z ~ 4, thereby determining its star-formation rate (S FR, via its 
FIR l uminosity, Lir, across rest-frame 8-1,000 fim, e.g. lKennicuttl 
Il998h and its characteristic dust temperature, T c \. Beyond the local 
Universe, these quantities have only rarely been well constrained 
until now. Alternatively, we can adopt a reasonable value for Td 
and then estimate the r edshifts of the SMG s seen in the vicinity 
of HzRG signposts (e.g. Stevens et al. 2003) us ing their rest-frame 
FIR c olours (e.g. lEales et alJl2003l;lGreve et alfeooijl : iPenner et all 
1201 lh . i.e. without resorting to conventional spectroscopy, thereby 
determining the likelihood that they inhabit a massive structure 
alongside the HzRG. 

In the next section we describe an extensive set of obser- 
vations, then present the reduced images, spectra and associated 
analysis in !j3] Sj4] contains our interpretation of those data and 
discussion of their implications. We finish with our conclusions 
in fj5] Throughout the paper we use a cosmology with Ho = 
71kms- 1 Mpc" 1 ,r2 m = 0.27 and ft A =0.73. 



z ~ 2 SMG field population ITacconi 

l2008l : lYounger" et alj |2008). From this we might conclude that the 
mechanism for the formation of the very massive galaxies in these 
over-dense regions may be funda mentally different to that of the 
>L* galaxies forming from SMGs dSmail et al.l2004 | ). although the 



earliest 12 CO J= 1-0 imaging of SM G s (e.g.lCarilli et alJl2oTTl : 



llvison et al.l201ll : lRiechers et alhoiTal R lHodge et al.l201lh using 
the Karl Jansky Very Large Array (JVLA) dPerlev et al.Eoi 1 ) sug- 
gests their cold gas may be more extended than the high- J 12 CO 
tracers employed previously. 

Valiant, early attempts to detect 12 CO to wards HzRGs with 
single-dish telescop es ended in failure (e.g. lEvans etlril 1 19961 : 
Ivan Oiik et alj[l997l) , but the stable spectral baselines afforded by 
interferometers eventually allowed the secur e detection of 12 CO 
J = 3 — 2 towards 53W002 at z = 2.39 dScoville etal]|l997l ; 
lAlloin, Barvainis & Guilloteaul |2000|) and 12 CO J = 4 - 3 to- 
wards 4C 60.07 at z = 3.79 and 6C 1909+72 at z = 3.53 
( Papadopoulos et al 1 l200d. hereafter POO), B3 2330+3927 at z = 
3.09 dDe Breuck et al.ll2003l hereafter DB03), TNJ0121+1320 
at z 3.52 dDe Breuck Neri & Omontl |2003|) and 4C41.17 at 
z = 3.80 dDe Breuck et al.l2005l) . D etections of the 12 CO J = 1-0 
line - that most sensitive to cold gas ( Papadopoulos & Ivisonl2002h 
- are particularly difficult towards radio galaxies: synchrotron 
emission often dominates at the frequency of the line, ~115 GHz, 
which is expected to be roughly an order of magnitude fainter 
tha n 12 CO J = 4 -3. The robust detection of 12 CO J = 1 — 
by lEmonts et alj d201 lh towards MRC 0152-209 is the excep- 
tion that proves the rule, albeit at a modest re dshift (z = 1.92) . 
More often, J = 1 — has proved elusive dlvison et al.1 1 19961 : 



The observations described hereafter targeted the HzRGs, 
6C 1909+72 and B3 J2330+3927. The former presents a clas- 
sic double-lobed (plus core) radio morphology, subtending 
~15 arcsec, and has a total 1.4-GHz flux density, 51.4GHz ~ 
259 mJy; the latter is brighter, Si.4GHz ~ 405 mJy, with what 
appears at first sight to be a similar morphology, if more com- 
pact (~2 arcsec), but which high-resolution radio imaging re- 
vealed to be an u nusually one-sided jet driven by a compact, 
flat-spectrum core dPerez-Torres & De Breuckl 12005b . Alongside 
PKS 1 138-262, 6C 1909+72 and B3 J2330+3927 are the m ost lu- 
minous of the 69 HzRGs explored by Sey mour et al.1 d2007l) in the 
rest-frame near-IR using Spitzer. 



2.1 JVLA observations 

Whilst the JVLA was in its most compact configuration during 
2011 September-November, we acquired amongst the first data 
taken in a new mode offering almost an order of magnitude more 
bandwidth than was previously possi ble, tuning to the 12 CO J = 
1-0 transition at 1 15.27120256 GHz dMorton & Noreaull 19941) for 
our target HzRGs. 

Short slots, usually 2-3 hr long, were scheduled dynamically 
to ensure excellent phase stability and transparency in the K and Ka 



1 Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA. 
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Table 1. 12 CO emission-line properties. 



Target 


R.A. 

(J2000) 


Dec. 
(J2000) 


Peak5„ 
(mjy) 


z 


FWHM 
(kms" 1 ) 


^CO 
(Jykms" 1 ) 


J = l-0: 
6C 1909+72 

B3 J2330+3927 (AGN host) 
JVLAJ233024.69+392708.6 (c) 


19:08:23.70±0.07 
23:30:24.69±0.02 


+72:20: 11. 63±0.36 
+39:27:08.58±0.16 


0.26 ±0.04 
3cr < 0.28 
0.21 ± 0.04 


3.5324 ± 0.0007 
3.0884 ± 0.0010 


570 ± 90 
720 ± 170 


0.222 ±0.049 

3cr < 0.074 
0.162 ±0.034 


J = 4-3: 
6C 1909+72 

B3 J2330+3927 (AGN host) 
JVLAJ233024.69+392708.6 (c) 


19:08:23.73±0.02 
23:30:24.85±0.02 
23:30:24.62±0.02 


+72:20:11.54±0.11 
+39:27: 12.04±0. 15 
+39:27:08.46±0.18 


3.3 ±0.3 
2.7 ±0.3 
1.9 ±0.4 


3.5324 ± 0.0006 
3.0934 ± 0.0005 
3.0901 ± 0.0006 


800 ± 90 
830 ± 100 
520 ± 110 


2.69 ±0.27 
3.29 ±0.51 
1.12 ±0.17 



atmospheric windows for 6C 1909+72 and B3 J2330+392, respec- 
tively. During these slots, data were recordec0 every 3 sec in each 
of 2 x 8 contiguous baseband (dual-polarisation) pairs, each base- 
band comprising 64 x 2-MHz channels for a total dual-polarisation 
bandwidth of 2,048 MHz, well over 20,000 km s _1 at the redshift 
of 6C 1909+72. Bright, compact calibration sources - flat-spectrum 
blazars lying within a few degrees of our target galaxies - were 
observed every few minutes to determine accurate complex gain 
solutions and bandpass corrections. 3C48 and 3C286 were also 
observed to set the absolute flux scale, and the pointing accuracy 
was checked locally every hour. 

For 6C 1909+72 the two sets of eight contiguous baseband 
pairs were themselves placed contiguously. We shifted the 12 CO 
J = 1-0 line (expected at 25.434 GHz for z = 3.5322) by 180 MHz 
from the centre of the available bandwidth to avoid the edge of a 
baseband, or one of the end basebands. 

For B3 J2330+3927, the 12 CO J = 1-0 line lies below 32 GHz 
in the Ka band. As such, only one set of eight contiguous baseband 
pairs could be deployed on the line (expected at 28.156 GHz for 
z = 3.094, but purposefully offset by 64MHz to avoid the edge of a 
baseband); the other set of eight basebands were tuned to 32.5 GHz. 

The data were red uced and imaged u sing ALVS following 
the recipes described bv llvison et aLl ( !201ll) . though with a number 
of significant changes: data were loaded using BDF2AIPS, avoid- 
ing any compression, and FRING was used to optimise the de- 
lays, in software, based on 1 min of data for 3C 48 or 3C 286. The 
basebands were knitted together using the NOIFS task, yielding uv 
datasets with either 1,024 or 512 x 2-MHz channels. Continuum 
images were made, using all of these data to obtain very deep maps 
of the synchrotron emission at 25.4 and 28.2 GHz (for 6C 1909+72 
and B3 J2330+3927, respectively), which was then cleaned and 
subtracted from the uv data (using UVSUB). The channels were 
then imaged in groups of four (yielding ~90kms _1 velocity res- 
olution) over ~3-arcmin-diameter fields, with natural weighting 
(ROBUST = 5), to form large cubes centred on the radio galaxies. 
The resulting spatial resolutions were 4.0 x 3.0 arcsec 2 at position 
angle (PA), 171°, for 6C 1909+72 and 3.4 x 2.8 arcsec 2 at PA, 72°, 
forB3J2330+3927. 



2.2 IRAM PdBI imaging 

Observations of 6C 1909+72 were made in 12 CO J = 4-3 us- 
ing PdBFs C and D configurations in the 3-mm atmospheric win- 
dow during 2007 July (with five antennas) and December (with six 
antennas) in a series of ~6-ll-hr tracks with good phase stabil- 
ity (0.6-1.0 arcsec seeing) and transparency (1.5-2.5 mm of water 
vapour). The local calibrator, 1928+738, was used to track ampli- 
tude and phase, with 1642+690 used for an independent check. Ab- 
solute fluxes were calibrated using MWC 349. Data cubes were 
made via ALVS using a natural weighting scheme, yielding a 
3.1 x 2.7 arcsec 2 beam with the major axis at a PA of 47°. 

The original PdBI data acquired bv lDB03l for B3 J2330+3927 
were re-imaged using ALVS using a natural weighting scheme, 
discarding those channels observed only in the compact configura- 
tion, yielding a 3.2 x 2.4 arcsec 2 beam with the major axis at a PA 
of 72°. 

2.3 Herschel imaging 

The Herschel data presented in this paper comprise a small part 
of the Hersc hel Radio Galaxy E volution Project (HeRGE) OT1 
programme dSevmour et al.l |2012|) . With PACS, we made near- 
orthogonal scans of our HzRG targets at nominal speed (i.e. with 
the telescope tracking at 20 arcsec s -1 ) using mini-map mode, 
spending a total of 288 sec on source, recording data simultane- 
ously at 100 and 160 /im, each scan comprising 12 x 4-arcmin scan 
legs separa ted by 4 arcse c. With SPIRE we approached the confu- 
sion limit dNguven et al.ll2010T) with just three repetitions (111 sec 
on source) of its 'small map' mode, obtaining data simultaneously 
at 250, 350 and 500 /im. The SPIRE images were produced using 
the standard pipeline; the PACS data were tackled with a variant of 
the pipeline developed bv llbar et alj ifeoiO). 

Flux densities for SPIRE were measured using SUSSEXTRAC- 
TOR JSmith et ap2012l) . as implemented in the Herschel Data Pro- 
cessing System. For PACS, we used apertures of radius 9 and 
13.8 arcsec at 100 and 160 /j,m, respectively, with appropriate aper- 
ture corrections, with errors determined by placing many random 
apertures in regions of the image with integration times similar to 
those of our targets. 



2 A bug in the Correlator Back End computer during this period resulted 
in only 1 sec of data being saved per integration, meaning that a significant 
fraction of the useful data were lost. 



3 RESULTS 
3.1 6C 1909+72 

The synchrotron continuum emission from 6C 1909+72 dominates 
the observed flux density at 25 GHz, being almost an order of mag- 
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Figure 1. Left: False colour image of 12 CO J =1-0 emission from 6C 1909+72, as measured by JVLA. N is up; E is left. Inset: 12 CO J = 4—3 as measured at 
IRAM PdBI, on the same spatial scale. Emission coming towards and going away from us relative to the systemic velocity is represented with the appropriate 
colours, revealing hints of a velocity field common to both 12 CO lines, though not globally. The much brighter 25-GHz synchrotron emission, against which 
we have struggled to discern the 12 CO J = 1 — emission, is represented with isophotal contours at —3, 3, 6, 12... X the local noise level. Right: 12 CO 
J = 1 — spectrum of 6C 1909+72, after subtraction of the dominant synchrotron emission component, smoothed with a 140-kms — 1 FWHM Gaussian. Inset: 
zoomed in on the line, with the new IRAM PdBI 12 CO J = 4— 3 spectrum ( i]2.2t shown in red, binned to 100 km s — 1 and scaled by 16 _1 X to be on the same 
Rayleigh-Jeans brightness temperature (Ti-,) scale as the 12 CO J = 1 — spectrum. Th e velocity scales correspond to z = 3.5322, approximately midway 
between the two components (z = 3.5203 and 3.5401) identified by Smith et alj d2010h . The red horizontal line shows the spectral region summed to create 
the 12 CO J = l-0 image shown alongside, and to determine Iqo- 
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Figure 2. Left: Greyscale K' image of the B3 J2330+3927 field (DB03) with isophotal contours of the 12 CO J = 1-0 emission in red (-3, 3, 3^, 6... X 
local noise level). N is up; E is left. The 12 CO J = 1 — is associated with compo nent c (circled, lower righ t), rather than the core of the AGN host galaxy 
(also circled, betrayed by its powerful synchrotron emission - cyan contours - Perez -Torres & De Breucla2005l) . 12 CO J = 4— 3 emission is shown in yellow 
- a radically different morphology from that presented by |PB03l . The brightest component is centred on the AGN while a fainter clump lies to the south-west, 
coincident with component c. The astrometric uncertainties here are ^,0.5 arcsec. Right: 12 CO J=l-0 spectrum of JVLA 233024.69+392708.6 (component 
c), near B3 J2330+3927. Inset: zoomed in on the line, with the IRAM PdBI 12 CO J = 4— 3 spectrum of the same galaxy ( >i2.2t shown in red, binned to 
53 km s — 1 and scaled by 16 — 1 X to be on the same Rayleigh-Jeans scale as the 12 CO J = 1—0 spectrum. The red horizontal line shows the spectral region 
summed to create the 12 CO J = 1 — image shown alongside, and to determine Ico- 
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Table 2. Measured flux densities and properties of HzRGs and brightest field SMGs. 



Target 


•SlOOfim 

(mJy) 


S'l60fim 

(mJy) 


J 250/Jm 

(mJy) 


(mJy) 


°500/jm 

(mJy) 


•S , 850nm 

(mJy) 


6C 1909+72 

SMMJ190827.5+721928 
B3 J2330+3927 
SMMJ233019. 1+392703 


17.0 ±3.0 
0.0 ±2.9 

11.3 ±2.9* 
5.4 ±3.3 


34.9 ±6.5 
-0.1 ±6.7 
27.0 ±6.2* 
4.2 ± 11.9 


57.2 ±2.7 
45.1 ± 2.6 
52.0 ±2.7 
31.6 ± 2.8 


69.6 ±2.8 
65.1 ± 2.7 

60.3 ± 2.7 

35.4 ±2.6 


63.4 ±3.3 

61.2 ±3.0 

57.3 ± 3.3 
21.1 ± 3.4 


34.9 ±3.0 
23.0 ± 2.5 
22.2 ± 2.7 
8.2 ± 1.9 



t Uncertainty does not include the contribution from confusion noise, nor any possible systematic offset in flux calibration. 
* Astrometry is consistent with a significant fraction of the PACS emission originating from JVLA 233024.69+392708.6 (Fig. [5). 



nitude brighter than the 12 CO line peak (Fig.[T| left). To assess the 
12 CO line properties we subtracted an accurate model of the contin- 
uum emission, exploiting the impressive bandwidth now available 
at JVLA. 

In the resulting continuum-free JVLA datacube, 12 CO J = 
1 — line emission is seen clearly, centred at z = 3.5324 ± 0.0007 
(Fig. right), closely matching that seen in our 12 CO J = 4 — 3 
spectrum from IRAM PdBI, z = 3.5324 ± 0.0006. The basic char- 
acteristics of the 12 CO line emission - peak and total flux densi- 
ties, line widths, etc. - are listed in Table Q] for both transitions. 
It is notable that we measure 7co4-3 — 2.69 ± 0.27Jykms~ , 
cf. the 1.62 ± 0.30 Jy km s _1 measured bv lPOOL who may also have 
under-estimated the line width due to a lack of available bandwidth. 

Both the 12 CO J = 1 — and J = 4 — 3 emission are par- 
tially resolved, even with a natural weighting scheme. Some veloc- 
ity structure can be discerned in Fig. [I] common to both transitions, 
but there does not appear to be a strong, coherent gradient. 

From the measured 12 CO properties, we find icoi-o = 
(1.19 ± 0.26) x 10 11 Kkms _1 pc 2 and the molecular gas mass, 
Mga S (H2 + He), is therefore likely around ~ lO n M0, 
naively adopting M^/L gq = 0.8 M Q (Kkms -1 pc 2 ) -1 
dDownes & SolomorJI 19981) . 



3.2 B3J2330+3927 

At the position of the radio galaxy core, as identified by 
IPerez-Torres & De Breuckl d2005l) and detected in 12 CO J = 4 - 3 
using PdBI by dDBolir we find synchrotron emission (S28.2GHZ 
5, 1 ± 0.02mJy centred at 23 h 30 m 24. s 866, +39°27'11.84" 
J2000) but there is no evidence of 12 CO J = 1 — emis- 
sion in a synchrotron-subtracted data cube. The (3-cr) limit we 
can set, assuming the same line width as the 12 CO J = 4 — 
3 emission (830 km s" 1 FW HM) and following Appendix A of 
ISeaquist. Ivison & Halll ( fl99l) is Inn < 0.074 Jykms -1 . The ba- 
sic characteristics of the 12 CO line emission are again listed in Ta- 
ble Q] for both transitions. 

Fig. |2] shows the Keck/NIRC AT-band imaging (from |PB03l) 
as a greyscale, superimposed with contours to show the emission 
in 8.4-GHz continuum (synchrotron), 12 CO J = 1 - and 12 CO 
J = 4 — 3. Approximately 4 arcsec to the SSW of the AGN (and 
~260kms -1 bluewards of its 12 CO J = 4-3 redshift) we do 
see emission in 12 CO J = 1 — 0. It is coincident with what lDB03l 
labelled component 'c' and tentatively suggested might lie at the 
same redshift as the radio galaxy. We show its spectra in 12 CO 
J = 1 - and J = 4 - 3 in Fig. |2] Our new 12 CO J = 1 - 
imaging and our re-analysis of the 12 CO J = 4 — 3 data make it 
clear that component c is in fact a gas-rich galaxy in the immediate 
vicinity of the radio-loud AGN, which is presumably interacting 
with (and thereby triggering activity in) the host galaxy of the AGN. 
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Figure 3. Three velocity slices in 12 CO J = 4- 3 towards B3 J2330+3927 
displayed as a false-colour RGB image, showing how the super-thermal 
emission centred near the radio-loud AGN (shown as black dots) is offset to 
the red from the 12 CO emission associated with the neighbouring, gas-rich 
galaxy, JVLA 233024.69+392708.6. Cyan, yellow and red contours (start- 
ing at 3 (j and spaced by \/2) show the 16 0-, 250- and 8 50-^tm continuum 
emission as seen by Herschel and SCUBA I Steve ns et alj2003l) , with no as- 
trometric tweaks applied, consistent with cold dust distributed throughout 
the system. N is up; E is left. 



Component c is detected in at least three IRAC bands and so has 
a significant stellar population. It is not obvious that it contributes 
significantly to the L yman a halo that surrounds B3J2330+3927 
dMatsuda etaT]|2009h . 

Fig.[3]reveals the dynamical structure of the 12 CO emission, 
with a blue-red gradient leading from component c to beyond the 
radio-loud AGN. The continuum emission at A b s = 160, 250 and 
850 fj,m (rest-frame 40, 60 and 200 /nm) seen by PACS, SPIRE and 
SCUBA are also shown. At A D b s = 24 /im, the emission is pre- 
dominantly from the AGN; at A b s = 100 and 160 ^m, it is pre- 
dominantly from component c; at A Q b s ^ 250 /im the emission is 
centred between the two. We have evidence, then, that dust at a va- 
riety of temperatures is distributed throughout the system, yet we 
lack the spatial resolution to disentangle these components reliably. 

The total flux in 12 CO J = 1 — for component c, 
^co = 0.162 ± 0.034 Jykms -1 . The excit ation situation for 
component c is c omparable to that of SMGs (Harri s~et al.ll2O10t 
llvison et all 1201 lh . with a brightness temperature (Tb) ratio, 
ico4- 3 /icoi-o = r-4-s/i-o = 0.43 ± 0.11 (cf. 0.41 ± 0.07 - 
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Figure 4. Colour-co lour plots for the H zRGs and their neighbouring 
SMGs, adapted from Amblard et al. (2010), together with a new plot de- 
signed to exploit our SCUBA 850-/im photometry. For galaxies/bands with- 
out robust detections in Table [2] upper limits are shown (calculated as 
S\ + 3<t\ ). for The c oloured backgrounds indicate the redshifts of model 
SEDs jAmblardetaI]|2010l) . The red crosses represent B3 J2330+3927 
(thick lines) and SMMJ233019. 1+392703, where the latter appears to lie 
at a lower redshift than the radio galaxy. The black crosses represent 
6C 1909+72 (thick lines) and SMM J190827.5+721928, which have re- 
markably similar colours, consistent with a shared redshift of ^,3. 



emission, r 4 _ 3 / 1 _ > 2.7 (3a), i.e. the emission is super-thermajf] 
(see 34~TV 

|PB03l assumed r4_3/i_o = 0.45 to extrapolate their 12 CO 
J = 4 — 3 measurement to J = 1 — and calculate a gas mass 
for B3 J2330+3927. We can see now that this will have resulted 
in an over-estimate of M gas for the AGN host galaxy. Our limit 
on Zcoi-o implies L' co < 2.9 x 10 10 Kkms -1 pc 2 and Af gas < 
2.3 x 1O 1O M where M gas /Leo = O.8M (Kkms -1 pc 2 )" 1 . 

For the gas-rich galaxy previously known as component c, 
which we re-name JVLA 233024.69+392708.6, we find L' co = 
(6.9 ± 1.5) x lO^Kkms^pc 2 and M gas ~ 5.5 x 10 10 M o . 

3.3 On the HzRGs and their environments 

The over-density of bright SMGs found around HzRGs, their faint- 
ness at optical/IR wavebands and the alignment of those SMGs 
with the jets from the central AGN were taken as tentative evi- 
dence that HzRGs are signposting massive structu res in which other 
galaxies are also un dergoing intense starbursts dlvison et alj|2000l : 
IStevens et al1l2003l) . 

Our deep Herschel imaging allows us to assess the likely red- 
shifts of th ese bright neighbouri ng SMGs via their FIR/submm 
colours (e.g. l Amblard et al.12 010. though such colours are only sen- 
sitive to (1 + z)/Td ~ see Blain 1999). We report their FIR/submm 
flux densities in Table|2] 

The upper panels of Fig. [4] show two of the colour-colour di- 
agnostic plots employed bv lAmblard et alj<201fjp to assess the red - 
shift and T d of galaxies detected by H-ATLAS dEales et alJ|20Tch . 
probing their colours across the rest-frame ~100-^im bump. The 
coloured backgrounds indicate the redshifts of model SEDs. The 
new, lower panel of Fig.|4]is based on the same model SEDs and is 
designed to exploit information from our SCUBA 850- /im imaging 
that is particularly relevant to galaxies at z > 2. 

The red crosses represent B3J2330+3927 and its brightest 
SMG neighbour, SMM J233019. 1+392703. The position of the lat- 
ter - particularly in the lower two plots - suggests strongly that it 
lies at a lower redshift than the radio galaxy, nea rer the peak of 
the ra dio-detected SMG population at z ~ 2.2 (e.g. Chapm an et al.l 
l2005h . 

The black crosses represent 6C 1909+72 and the very bright, 
nearby SMG, SMM J190827. 5+721928. Their rest-frame FIR 
colours are remarkably similar in the lower two colour-colour plots. 
The two galaxies differ in S160/S250, ie- the colours on the Wien 
side of their SEDs. The limit for the bright companion is consistent 
with the known redshift of the radio galaxy, but the colour of the ra- 
dio galaxy is bluer than we might expect. On balance, the evidence 
is consistent with the proposition that these two dusty starbursts 
share the same node or filamenjf] of the cosmic web. 

With the Av > 10, 000 km s _1 of velocity coverage available 
to us via the WIDAR correlator at JVLA, we should be sensitive 
to 12 CO J = 1 — emission from bright SMGs in our target fields, 
should they have a similar ratio of Sco / Sgso/jm an ^ ^ ~ as seems 
likely, at least for 6C 1909+72 and SMM J190827.5+721928 - they 
lie in the same cosmic structure as the central HzRG. 



Both well et alj|2012h . This differs strikingly from that of the radio 
galaxy core, where the Tb ratio is entirely inconsistent with thermal 



3 Strictly, super-thermal excitation exists only where local thermodynamic 
equilibrium has been violated by a population inversion. Here, we use it as 
short-hand for r > 1. 

4 A typical filamen t at z ~ 3.5 spans of order 30 h~ 1 Mpc (e.g. 
ISprineel et afll2005t) . SMM J190827. 5+721928 is separated by 45arcsec 
(335 kpc) in the plane of the sky from 6C 1909+72. It would need to lie 
within Sz 0.04 to inhabit the same sheet or filament. 
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Figure 5. 850-^m continuum emission (\/2-spaced contours, starting at 3 a, IStevens et alj|2003l) for the fields surrounding B3J2330+3927 (left) and 
6C 1909+72 (right), superimposed on false-colour images made using the Herschel 160-, 250- and 350-/^m imaging, described in i]2.3l convolved to the 
350-/im spatial resolution. N is up; E is left. The SMGs identified bv lStevens et"ai] J2003I) are circled. Unsurprisingly, these are revealed here as the reddest 
objects in each field. The red structure apparent in the 6C 1909+72 field is co-aligned with both the radio jets (recall Fig.[T} and a prominent extension seen 
previously at 850 fim. A less prominent extension is seen to the WNW of B3 J2330+3927. 



We do not detect CO J = 1 — emission from 
SMMJ233019. 1+392703 in the B3J2330+3927 field. However, 
this SMG is less than half as bright as the nearby HzRG its emis- 
sion is attenuated severely (82 per cent) by the JVLA's primary 
beam at 28.2 GHz, so this does not rule out the possibility that it 
shares the same structure at z ~ 3.09. We set a 3-cr upper limit of 
/coi-o < 0.41Jykms _1 where we have adopted a line width of 
800 km s -1 , typical for SMGs (e.g. lGreve et~ai1l2005b . 

SMMJ190827.5+721928, in the 6C 1909+72 field, is suffi- 
ciently bright that we would expect to detect its 12 CO J = 1 — 
emission, given our Av ~ 20,000kms _1 [3.39 < z < 3.75] of 
velocity [redshift] coverage in this field, if it were to have a similar 
ratio of Scoi-o/SssOfim an d be part of the same structure as the 
nearby HzRG. We searched the 12 CO J = 1 — velocity cube at 
the position of the SMG, at spectral resolutions ranging from 95- 
380 km s~ 1 . Despite being amongst the brightest known SMGs, its 
position is ill-defined since it lacks a convincing counterpart in the 
available Spitzer 3.6-24-/im imaging. The most convincing peak 
lies at !/ b s = 25.124 ± 0.006 GHz, which would correspond to 
icoi-o = 0.110±0.033 Jykms- 1 a!tz = 3.588±0.001 (after pri- 
mary beam correction). The best-fit line width (430 ± 120 km s^ 1 ) 
would be narrower than the majority of SMG 12 CO J = 3-2 lines 
dGreve et"al . 2005; Bot hwell et al.l2012h and we do not view this as 
a robust detection, but it would be an obvious place to begin search- 
ing for 12 CO J = 4—3 with IRAM's new WideX correlator. Several 
other 2— 3-cr peaks are also present, but non e with a line width even 
remotely commensurate with an SMG (e.g. lGreve et al .120051) . The 
1-f Ico level for a line of width 800 km s~ x is 0.06 Jy km s _1 . We 
are forced to conclude that the presence of two unusually bright 
SMGs in this field may be due to the chance alignment of distant 
starbursts. 

The lack of a convincing detection of 12 CO towards the 
brightest SMG in the vicinity of 6C 1909+72 casts doubt on the 



idea that it co-habits a proto-c luster environment with the ra- 
dio galaxy dStevens et al.l [2003). However, our Herschel imag- 
ing does hint at the presence of a dust-rich strand of the cosmic 
web centred on 6C 1909+72. Fig. [5] reveals a large (~500kpc) 
red structure, the majority of which is co-aligned with a promi- 
nent extension (~150kpc, somewhat larger than the ~50-kpc 
scale of th e northern radio jet - recall Fig. Q] seen previously 
at 850 dStevensetai]|2003l) . We might imagine this as a fil- 
ament of gas that is being accreted slowly onto the AGN host 
galaxy, but in our view it is no less likely that this structure is 
due to the outflow of metal-rich material from the radio galaxy 
host, heated by a network of unresolved LIRGs. A dusty out- 
flow could be driven by radiation p ressure from a strong nuclear 
starburst and/or AGN on dust (e.g. [P rochaska & Hennaw i] BOO 1 
Lehnert et al. 2011; Faucher-Giguere & Ouataert 2012; Roth et ... 
120121; IWagner. Bicknell & Umemural 1201 2I) with material trans 
ported well beyond the AGN host galaxy, as witnessed in the 
near-IR on smaller scal es for radio g alaxies at z ~ 2 — 
3 jNesvadba et"ai] 120081; iHarrison et al] |2012h. for M87 in th e 
Virgo cluster dSimionescu et al.ll2008l 120091 ; IWerner et al.l l2010h . 
for the brightest cluster galaxy (BCG) in the Hydra A cluster 



( Kirkpatric k et al.ll2009l; 



erally ( Kirkpatrick et al 



iGitti et alj201 ll) and for BCGs more gen- 
120111; lO' Sullivan et al.ll20li1) . We note 
that the presence of ~3,000-kms _1 abs orpti o n featu rejfl against 
an unpolarised optical continuum led IPevI dl997t) to classify 
6C 1909+72 as a broad-absorption-line radio galaxy, and is con- 
sistent with a powerful outflow. Of course, we must also bear in 
mind the well-rehearsed argument that radio galaxies are most eas- 



5 Fora~3,000-kms _1 outflow, with a dynamical timescale of a few xlO 8 
yr, a single epoch of starburst and/or AGN activity would be sufficient to 
generate a structure on the scale seen here. 
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Figure 6. Radio-through-optical SEDs of the radio galaxies, B3 J2330+3927 (left) and 6C 1909+72 (right). Note that all quantities are observed rather than 
rest-frame. Their entire 3 fim to 74 MHz SEDs can be fitted adequately (see 94. It using only thermal dust and synchrotron emission components, whose 
characteristics are listed in Table [3] Our new measurements from PACS and SPIRE - which probe the SED peaks, constraining Lir and Td - are shown in 
red. Lacking the spatial resolution to account for contamination of flux densities for B3 J2330+3927 by component c (JVLA 233024.69+392708.6), we fit to 
the summed flux densities of both galaxies. 



ily d iscovered when their jet s are encountering a working surface 
(e.g. lBarthel & Arnaudll 19961) . 



4 DISCUSSION 

4.1 On the HzRG SEDs and T b ratios 

Fig. [6] shows the observed radio-through-optical SEDs of 
6C 1909+72 and B3 J2330+3927. For both galaxies, the photom- 
etry between 3 /im and 74 MHz can be described adequately with 
a model comprising on l y syn chrotron and thermal dust emission, 
following IKovacs et alj ^20 id ). Table[3]lists the results of a simul- 
taneous fit to the following free parameters: the dominant cold dust 
temperature, Ta, and the power-law index, 7, of a dust tempera- 
ture distribution, dM^/dT^ oc T^ 7 , designed to offer a physi- 
cally motivat ed treatment o f the Wien side of the thermal emission 
spectrum (in IKovacs et all T c = Td, the low-temperature cutoff 
of the distribution), the dust mass, Ma (for a characteristic pho- 
ton cross-section to mass ratio , Kssopm = 0.15 m 2 kg -1 , from 
iDunne. Eales & E dmunds 2003, where the frequency dependence 
of the dust emissivity, /3, was fixed to +1.5) and the synchrotron 
power-law in dex, a. The resulti ng measurements of Lir and giR 
(as defined bv lHelou et alj 19851 , but where Sir is measured across 
Arcst = 8— 1, 000 /J,m, as is Lir) are also listed. 

Both radio galaxies have warmer dust temperatures, T RzRG ~ 
45 K, than those seen for similarly luminous, dusty starbursts at 
z ~ 3 (T^ MG ~ 35 K when calculat ed using a power-law tem- 
perature distribution, as we have here - iMagnelli et al.ll2012l) . De- 
termining the contribution of the AGN to Lir is a long-standing 
problem that we cannot hope to solve without dramatically im- 
proved spatial resolution. The relatively warm dust does imply sig- 
nificant AGN contributions to Lir for these two radio galaxies, 
though even this cannot be asserted with certainty since the bright- 
est 100-160-Atm emission in the B3J2330+3927 system ( fl3~2l l 
is centred near the gas-rich companion (which may contain its 
own AGN, if it i s anything like the interacting system seen to- 
wards 4C 60.07 - Ilvisonet alll2008l) . On the other hand, the de- 



Table 3. T b ratios and parameters from spectral energy distributions. 



Parameter 


6C 1909+72 


B3 J2330+3927 


r 4-3/l-0 


0.76 ± 0.18 


3cr > 2.7 




1.14 


1.00 


7d(K) 


45.7 ±1.3 


41.4± 1.1 


logA/ d (M Q ) 


9.37 ± 0.04 


9.26 ±0.17 


7 


5.51 ± 0.05 


5.54 ±0.04 


a 


-1.22 ±0.04 


-0.95 ± 0.07 


<?IR 


-1.02 ±0.04 


-0.40 ±0.11 


log Lir (Lq) 


13.69 ±0.08 


13.52 ± 0.07 



tection of dust and molecular gas towards the radio galaxy hosts 
makes it likely that significant levels of star formation are tak- 
ing pl ace in both. Gi ven the strong dependence of Lir on Td 
(e.g. lEales et aDl2000h . we might expect that the fraction of Lir 
due to the AGN, Lfg N /L m « 1 - 7?(T d HzRG /Td SMG )- 6 , where 
77 w 0.5 is the frac tional AGN contri bution to Lir estimated for 
typical SMGs (e.g. iFraver etal.lll998l) . This suggests a ^90 -per- 
cent AGN-related contribution to Lir for our HzRGs, which still 
leaves room for vigorous sta rburst activity in the radio galaxy hosts, 
with SFR ~ 500 M Q yr" 1 dKennicuttll 19981) . 

The super-thermal Tb ratio, r4_3/i_n > 2.7(3(7), observed 
towards the AGN host galaxy, B3 J2330+3927, reveals the pres- 
ence of highly excited molecular gas, given the excitation re- 
quirements of 12 CO J = 4 — 3 (E4-s/k ~ 55 K and n C iit ~ 
1.9 x 10 4 cm" 3 ). Even for 6C 1909+72, the T b ratio, r 4 _ 3/1 _ = 
0.76 ± 0.18, is significantly higher than the average value seen 
for SMGs by iBothwell et at] d2012l) . IRAM PdBI and VLA ob- 
servations of 4C 4 1.1 7 at z = 3.80 revealed a similar st ory, with 
r 4-3/i-o likely super-thermal (Papadopoul os et al .120051) . In com- 
parison, the same line ratio averaged over the entire M ilky Way is 
r 4 _ 3 / 1 _o ~ 0.1-0.2 jFixsen. Bennett & MatheJ 19991) , indicative 
of quiescent gas in the interstellar medium (ISM). 

In contrast to B3 12330+3927, the T b ratio that we find for 
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its companion galaxy, JVLA 233024.69+392708.6 (r 4 _ 3/1 _ = 
0.43 ± 0.11) is consistent with the average value seen for SMGs 
and with star-forming regions within nearby starburst galaxies 
such as M82 and NGC253 , which have va lues in the range 
~0.5-0.8 dGtisten et al]|l993l ; iMao et alj|200d) . The T b ratio in 
JVLA 233024.69+392708.6 is thus consistent with gas heated by 
UV radiation from star formation, i.e. less extreme ISM conditions 
than those of the radio galaxy. 

Tb ratios as high as that observed towards B3J2330+3927 
are expected in environments where the gas energetics are dom- 
inated by X-rays. These can penetrate large columns of gas 
and maintain a high kinetic gas temperature, Tk, much deeper 
into molecular clouds tha n far-u lt raviolet (FUV) photons (e.g. 
Meijerink, Spaans & Israel 2006|; ISchleicher. Spaans & Klessea 
2010; van der Werf et alfcOlOh . with this gas thermally decoupled 



from the cooler dust. Powerful radio-loud AGN are often found 
to have high X-ray luminosities, capable of powering extended X- 
ray-dominated regions (XDRs); while B3J2330+3927 has no X- 
ray measurements, this is one possible explanation for the observed 

r 4 _ 3 /i_ . 

Another possible explanatio n is shock-excited gas (e.g. 
I Flower & Pineau Des Foret sl2010h caused by its radio jet ramming 
into the molecular environment, bearing in mind that the mechan- 
ical power of a radio jet can be orders of magnitude larger than 
that of its synchrotron emission. An example of this is observed 
in the powerful radio galaxy, 3C 293, which has Milky Way-like 
levels of star formation and a very low X-ray luminosity (see also 
iMatsushita et alj 2004, regarding the central region of M51). Its 
super-thermal 12 CO (as measured in 12 CO J= 1 — and J = 4— 3) 
has been attributed to turbulent heating from the dissipation of 
shocks caused b y the interaction of its powerful radio jet with 
its molecular gas jP apadopoulo s et al . 2008; lNesvadba et alfeoid : 
iGuillardetalfeoij) " - an inte raction that seemin gly also drives a 



~l,400-km s" 1 outflow in H I dMorganti et alj2003h . Such mecha- 
nisms are often proposed to provide the feedback needed to recon- 
cile models of galaxy formation with a variety of observations (e.g. 
Granato et aHl2004lM~cNamara & Nulsenll2012llFabianll2oT2l and 
references therein). 

A large fraction of the radio emission from B3J2330+3927 
comes from its core (~50 per cent at 8.4 GHz) and it displays an 
unusually one-sided jet: 'one of the m ost asymmetric radio struc- 
tures ever reported for a type II AGN' I Per ez-Torres & De Breuckl 
20051). This may be due to relativistic Doppler beaming (e.g. 

t ! 



Riley & Warnerl 1 1990) or because of the launch mechanism (e.g. 
Chagelishvili, Bodo & Trussoni 1996) but if a survey of such 
galaxies were to find relatively high Tb ratios relative to galaxies 
hosting two-sided jets then it would suggest a more prosaic reason: 
a blockage that causes an unusually high fraction of the mechanical 
energy from AGN-driven outflow to be deposited into the molecu- 
lar gas. 

According to IWagner & Bicknelll J20T l]), the kinetic energy 
and momentum of a jet can be transfered to dense gas with a high 
efficiency (10-70 per cent), causing turbulent heating and shock 
excitation of the molecular gas and - in extreme cases - driving 
gas outflows and inhibiting star formation. The kinetic en ergy of 
the rad io jet in B3 J2330+3927 can be estimated following Punsly 
d2005h and IPunslv et alj d2008h . assuming Sisim hz = 0.90 Jy 
(from S 3 65mhz = 0.394 Jy and a?t^ z = -0.93lDouglas et al.1 



1994 ICohen et al.ll2007l though see also iBirzan et alj|2008l and 



Cavagnolo et al.1 12010). The jet's kinetic luminosity, Q, is then 
~ 2 x 10 46 erg s" 1 , roughly two orders of magnitude larger than the 
12 CO J = 4-3 luminosity, L' C04 _ S = 3.4x 10 44 ergs" 1 . This sim- 



plistic approach assumes that the 151-MHz emission comprises op- 
tically thin radiation from the lobes, whereas a substantial fraction 
may come from a Doppler-boosted jet. Nevertheless, over the ~10- 
Myr lifetime of the radio source we can expect a few 10 60 erg to be 
deposited into the multi-phase ISM (enough to quenc h a cooling 
flow for ~ 10 9 yr, as argued bv lMcNamara et al.12 005). The turbu- 
lent kinetic energy of the molecular gas can be calculated through 
El1n b = l-MHa^Ha 2 , where cth 2 = FWHM/2\/21n2 is the ve- 
locity dispersion of the gas. From $3?2\ M H2 < 2.3 x 10 10 M 
so if (Tcoi-o ~ o"co4-3 = 350 km s" 1 and is solely due to the 
velocity dispersion of the gas (i.e. no components due to rotation 
or bulk outfl ow), then Ejf* h < 1 59 erg. We can thus conclude, 
as did (as did Nesvadb a et alj2008l for an analagous situation) that 
the radio jet carries sufficient energy to explain the observed CO 
characteristics. The jet can influence the gas properties in precisely 
the way required to stop star formation, though questions clearly 
remain about the ultimate fate of the gas and deeper observations 
are require d to search for an y high-velocity molecular component 
(see, e.g., |Polletta et alj|201 lh . 

Taking another approach, if we adopt a fiducial jet ad- 
vance speed, Vj c t ~ 0.1c, then the gas is heated to Thock ~ 
3/16 n mn/k vf et w 10 10 K, where fj, is the reduced mass and k is 
Boltzmann's constant. For a large fraction of the gas that has been 
shock-heated by the jet, th e density will be low (10~ 2 -10~ 3 cm" 3 , 
Wagner & Bicknell 2011). This over-pressurised gas will expand 
dBegelman & Cioffill 19890 . running into the denser, ambient ISM, 
which will be engulfed and de stroyed by the resulting shock (e.g. 
iKlein. McKee & Colell Jl 1994 ). The velocity of the shock that is 
driven into the dense clouds, « s , c i ou d = Uwind(Pwmd/ / Pcioud) 1/ ' 2 , 
where p„ ln d is the density of the hot plasma created by the 
outward-moving radio jet and p c i ou d is the density of the engulfed 
ISM. For typical molecular clou ds, p c i ou d ^10 4 cm 3 , so we ar- 
rive at v B , doud ~ 20^0 km s" 1 dBegelman & Cioffill 19891) . Such 
molecular shocks would be significant emitters in high- J CO lines 
dFlower & Pineau Des Foretsll2010l) . 

At what rate must the gas must be shock-heated to main- 
tain this CO emission? For a C-type shock wave with « S)C ioud ~ 
40kms~ 1 (assuming a lower velocity or a denser ambient ISM 
will make little difference, and a J-type shock will produce a line 
luminosity almost two orders of magnitude fainter), we find that 
to explain L'cqa-^ requires a mas s shock rate of ~ 30 Mq yr -1 
(Fl ower & Pineau Des Foretsll2010l) . Over the age of the jet, only a 
few x 10 9 Mq of gas needs to be shock heated, perhaps much less 
if this is a special, short-lived phase due to the initial passage of 
the radio jet. The gas may return to a higher density following the 
passage of the shock dCooper et al1l2008h . perhaps re-forming H2 
dGuillardetai] |2009). so the amount of shock-heated gas required 
to explain the observed Tb ratios is not restrictive. 

Whatever the excitation mechanism, the 12 CO optical depth 
must be moderate to low (ti_o ;S 1) since in the optically thick 
regime the Tb of a transition is equal to the excitation temperature, 
i.e. Tb ~ T ex (w ith the latter being the sam e for all levels when 
thermalised - see Papadopoulos et al1l201ll) . In the optically thin 
case, y*- 3 /y£-o x T 4_3/ti_o and since tj+i.j oc ( J + 1) (at 
least up to J ~ 7), we find T^~° < T 4 " 3 (true for thermalised gas, 
but also for sub-thermal gas if the effect due to optical depth is suf- 
ficiently large). Optically thin 12 CO emission can be achieved if the 
medium is highly turbulent, as is the case for the Galactic Centre - a 
local region where the 12 CO ladder is observed to be super-thermal. 
The super-thermal Tb ratio observed towards B3J2330+3927 is 
thus likely an indication of mechanical energy deposited into the 
ISM from the radio jet and/or the ongoing starburst, inducing a de- 
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gree of turbulence sufficient to make the 12 CO lines optically thin, 
allied with a high density n C rit) and a relatively high kinetic 
temperature; moreover, this must happen globally, rather than in 
some sm all fraction of the gas, as seen in the nucleus of M 82, for 
example dKnapp et alJll980| - |Weifi et alj|200ll) . 



4.2 The merger fraction amongst starbursting HzRGs 

The degree to which various populations of star-forming galaxies 
are driven by violent mergers of gas-rich galaxies or by the accre- 
tion of cold gas streams - sometimes smoo th, sometimes clumpy - 
via da rk matter filaments into spiral disks (IDekel. Sari & Ceverinol 
2009) has been a topic of intense debate in recent years. Even for 
SMGs - the most extreme star-formation events in the Universe, 
where the evidence for merger-driven ac tivity was thought by many 
to be overwhelming (Engel et al. 2010) - the seeds of doubt have 
been sown. 

What of radio galaxies? We cannot generalise to the entire 
population, even though that population is small in number: the 
time-consuming interferometric observations of molecular gas that 
are capable of laying bare the processes leading to star formation 
have generally been obtained towards systems from which rest- 
frame FIR emission has already been detected, i.e. we are limited to 
systems pre-selected to be undergoing intense starbursts. Even so, 
it is interesting to add our findings to those studies of z > 2 radio 
galaxies in the literature and explore the overall statistics. Taking 
the eight systems listed in fJT]in order of decreasing redshift: 

• The interfe rometric 12 CO J = 4 - 3 study of 4C41.17, at 
2 = 3.80, by |Pe Breuck et alj d2005l) revealed activity driven 
by a merger between two massive, gas-rich galaxies separated by 
~2 arcsec (~ 13 kpc) and ~400 km s - . 

• Hi gh-resolution interfe rometric imaging of 4C 60.07 at z = 
3.79 by llvisoiTet al. (2008) revealed two components of roughly 
equal integrated flux, separated by ~30kpc - evidence of an early- 
stage merger between the host galaxy of an actively-fueled black 
hole (the HzRG), and a gas-rich starburst/AGN, caught prior to its 
eventual equilibrium state. T his separation is typ ical of SMGs with 
multiple radio identifications dlvison etal .120071) and for the intense 
burst of star formation near first pa ssage in merger simulations (e.g. 
ISpringel. Pi Matteo & Her nquist 200J). 

• ISmith et alj ( 120101) use rest-frame UV spectroscopy to show 
that 6C 1909+72, at z — 3.53, comprises two galaxies separated 
by ~l,300kms _1 along the line of sight, with the radio-loud AGN 
hosted by the more distant component. 

• TNJ0121+1320 at z = 3.52 was detected in 12 CO J = 4- 3 
bv lDe Breuck et alj d2003h . There were no obvious signs of a gas- 
rich companion, but the available sensitivity, velocity coverage and 
spatial resolution were poor compared to recent studies. 

• We have shown here that the immediate environment of 
B3 J2330+3927 at z = 3.09 contains at least two gas-rich com- 
pon ents, again separated by ~30kpc. 

• iNesvadba et alj J2009) presented interferometric 12 CO imag- 
ing of the z — 2.58 radio galaxy, TXS 0828+193, finding two large 
gas reservoirs separated by 720 kins' 1 along a shared line of sight 
roughly 80 kpc SW of the AGN host galaxy, which was itself de- 
tected only in synchrotron at 96.6 GHz. 

• 53W002 at z = 2.39 was observ ed and detected in 12 CO 
dYamada et al.ll995MScovilre et all 19971) without knowledge of its 
FIR luminosity, though it was known to sit a mongst a clump of 
Ly-a emitting galaxies ( Pascare lle et al.|[l996l) through which the 
radio galaxy is suspected to be assembling a large fraction of its 



eventual stellar mass through mergers (e.g. lMotoharaetai1l200lh . 
At least one of the se - 330 kpc and £300 km s" 1 from 53W002 
- is FIR-luminous i Smail et alj|2003l) . Based on their 12 CO imag- 
ing, [Scovilleetal] 1 1997 ) concluded 'it is cle ar that the emissio n 



is resolved', but this was not corroborated by lAlloin et al.l (|2000|) . 
However, with a ~6-arcsec synthesised beam the latter data would 
be insensitive to close companions. 

• Finally, Hubble Space Telescope near-IR imaging of 
MRC 0152— 209 at z = 1.92 r evealed a morphology consistent 
with an advanced-stage merger ( iPentericci et alJuOOll) where we 
would expect the gas to have been driven into a compact configu- 
ration, well within the 7-10- arcsec FWHM synthe sised beam of the 
available 12 CO observations dEmonts et alj|201 lh . 

Overall, we have six mergers or close galaxy pairs, one am- 
biguous case, and one case where no merger activity can be dis- 
cerned (albeit with relatively poor sensitivity and resolution avail- 
able). The statistics suggest, then, that violent interactions are as 
ubiquitous amongst starb ursting high-red shift radio galaxies as 
they are amongst SMGs dEngel et aU2010l) . 



5 CONCLUSIONS 

Combining our new interferometric 12 CO data for two distant radio 
galaxies with those for six more, reported elsewhere, we find that 
activity in starbursting radio-loud AGN at high redshift is driven by 
the merger or interaction of two or more systems in which signifi- 
cant masses of molecular gas and stars have already formed, rather 
than the steadyff] accretion of cold gas from the cosmic web. 

We introduce a new colour-colour diagnostic plot to constrain 
the redshifts of several distant, dusty galaxies found previously in 
our target fields. We conclude that the SMG south of 6C 1909+72 
likely shares the same node or filament, but we fail to detect this 
FIR-luminous galaxy in 12 CO, despite ~ 20,000 km s _1 of velocity 
coverage. 

We introduce a new colour-colour diagnostic plot, exploiting 
the power of 350-850-/im photometry to constrain the redshifts of 
several distant, dusty galaxies found previously in our target fields. 
We conclude that the bright SMG near 6C 1909+72 likely shares 
the same node or filament of the cosmic web at 2 ~ 3.5 as the 
signpost AGN. However, we fail to detect this SMG in 12 CO J = 
1 — 0, despite our ~ 20,000 kms^ 1 of velocity coverage. 

Also in the 6C 1909+72 field, we find an unusually large, red 
dust feature, aligned with the radio jet. We suggest that metal-rich 
material may have been dispersed on ^100-kpc scales by a col- 
limated outflow, reminiscent of the jet-oriented metal enrichment 
seen in X-ray observations of local cluster environments. 

We find that 12 CO brightness temperature ratios in the host 
galaxies of radio-loud AGN are significantly higher than those 
seen in similarly intense starbursts where AGN activity is less pro- 
nounced. In our most extreme example, the super-thermal Tb ratio 
suggests that significant energy is being deposited rapidly into the 
molecular gas via X-rays and/or mechanical ('quasar-mode') feed- 
back from the AGN, leading to a high degree of turbulence globally 
and a low optical depth in 12 CO - feedback that will lead to the ces- 
sation of star formation on a timescale commensurate with that of 
the jet activity, £ 10 Myr. 



e Often denoted 'secular, which in astronomy indicates that a process 
occurs continuously, as opposed to a discrete or periodic event such as a 
merger. 
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